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It is thoroughly documented and established that noise 
can play a constructive role in neuronal systems. Noise 
correlations, which have been observed in the brain, are 
thereby usually assumed to be ignorable. However, it has 
been shown that noise correlations cannot be avoided and 
may indeed play a vital role in neuronal dynamics, mainly 
because they affect the amount of information transmit- 
ted across the cortex as well as the computational strate- 
gies of neuronal networks. In the context of transmission 
of neuronal information, there is an ongoing debate about 
whether a cortical neuron is sensitive to the mean firing 
rate of presynaptic neurons and their correlations or not. 
Regardless of this, the mean firing rate and the correla- 
tions among neuronal groups are two important factors 
determining the transmission of information across the 
network. Here we elaborate on the effects of correlated 
Gaussian noise (noise correlation) on the mean firing rate 
and the correlations among neurons of an electrically cou- 
pled neuronal network. We find that the noise correlation 
has little effect on the variations of the mean firing rate 
with respect to the noise intensity. Variations of the pop- 
ulation coherence with respect to the noise intensity, how- 
ever, are strongly influenced by the ratio between local and 
global Gaussian noise in the overall noise intensity, which 
can be tuned by the noise correlation. Our results indicate 
that noise correlation may have a significant impact on the 
response of postsynaptic neurons if these are sensitive to 
correlated neuronal activities. 



I. INTRODUCTION 



Neurons are usually subject to random fluctuations on dif- 
ferent scales, ranging from channel noise created by random 
ion flow across the plasma membrane to synaptic noise cre- 
ated by the activity of other neurons. In experiments, it has 
been shown that noise has constructive effects on neuronal 
dynamics. For example, William and Durand 1 showed that an 
appropriate noise intensity can improve the detection of sub- 
threshold signals in a resonant manner. Higgs et al. 2 found 
that synaptic noise increases the gain in many pyramidal neu- 
rons with large sAHP (slow after hyperpolarization). Jacob- 
son et al, 3 found that channel noise contributes significantly 
to membrane voltage fluctuations at the subthreshold voltage 
range. In theoretical and computational studies, the construc- 
tive role of noise in neuronal systems has been reported as 
well. For example, it was shown that noise is able to evoke 
coherence and stochastic resonance in single neurons^— as 
well as in one-dimensional 910 and two-dimensional^— neu- 
ronal networks. Related to the present study, in the sense that 
correlated noisy inputs have been considered, are the two pa- 
pers by Kreuz et al. 7,8 where it has been shown that an inter- 
mediate noise intensity can evoke the most coherent temporal 
output of a single neuron irrespective of the noise correlation 
length. Double coherence resonance in terms of an optimal 
combination of noise intensity and correlation was reported 
as well. Notably, the effect of auto- and cross-correlations 
of input spikes on the response of spiking neurons has also 
been studied extensively. 15 Additionally, noise can also induce 
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and/or enhance the synchronization in neuronal systemsjiiais 
Refs. O and [l8| are two comprehensive review papers, rec- 
ommended to the readers who are interested in the research of 
noise effects on nonlinear systems in general, including neu- 
ronal dynamics. 

In the cortex, a single interneuron can be connected with 
tens of thousands of local circuit interneurons. Thus, activ- 
ity of this single interneuron will provide correlated input 
to many neurons in the local circuit. 19 Furthermore, corre- 
lated inputs might be stimulus driven and usually are random. 
These random correlated inputs can be described by means 
of correlated noise. In the past, noise correlations have, apart 
from a few exceptions, 7 ' 8 usually been ignored. But some re- 
searchers have found that noise correlations can have many 
different effects on the neuronal population — the amount 
of information encoded, the computational strategies of net- 
works of neurons, etc., as reviewed in Ref. |20j This means 
that noise correlations should not be neglected in neuronal 
systems. Adding to this conclusion is also the fact that the 
conceptually related coupling via noise in one-dimensional 21 
and two-dimensional 22 systems has been found to induce syn- 
chronization. Meanwhile, the question of how a postsynap- 
tic neuron is affected by the presynaptic neuronal population 
is still not fully understood. There is an ongoing debate on 
whether a cortical neuron is driven mainly by the mean fir- 
ing rate of presynaptic neurons or by correlated firing activi- 
ties.— _ — No matter what the outcome of the debate, the mean 
firing rate and correlations of neuronal groups are two im- 
portant factors in investigating transmission of neuronal in- 
formation. Therefore, we will investigate the effects of cor- 
related Gaussian noise (noise correlation) on the mean firing 
rate and correlations of a neuronal population in this paper. 
The obtained results may have important implications for un- 
derstanding the transmission of neuronal information. 

The paper is organized as follows. Equations governing the 
two-dimensional neuronal network are presented in the next 
section. Measures used for quantifying the observed neuronal 
dynamics are introduced in section [HI] while the results due 
to correlated Gaussian noise are presented in section [IV] Fi- 
nally, the summary is given in section [V] We also provide an 
algorithmic description of noise generation in the Appendix. 



where u n (i, j) is the membrane potential of neuron and 
v n (i, j) is the corresponding ion concentration at the discrete 
time n. The system parameters are a, /3 and 7, whereby the 
latter two determine the time scale associated with the dynam- 
ics of the slow variable v n (i,j) and a is the main bifurca- 
tion parameter. If not stated otherwise, we use a = 1.99 
and (3 = 7 = 0.001, for which each neuron is governed by 
a single excitable steady state (u*,v*) = (—1, —1 — a/2). 
Each neuron is coupled electrically with its four nearest neigh- 
bors with periodic boundary conditions given by u(0,j) — 
u(N,j),u(N + l,j) = u(l,j),«(i,0) = u(i,N),u(i,N + 
1) = u(i, 1). Finally, D is the coupling strength between the 
neurons on the 128 x 128 spatial grid. 

The correlated Gaussian noise rj n (i, j) is expressed as 

Vn(i,j) = V~Re n + Vl - R£n(i,j), (1) 

where e„ is the Gaussian white noise and common to all units, 
i.e. global noise, with the properties: 

j (e„) = 0, (2) 

and £ n (i, j) is the local Gaussian noise, which is uncorrected 
from site to site. £ n is taken as Gaussian white noise with 
the properties: 

and Gaussian colored noise with the properties: 

f (f»(i,j)>=0, ; 

I (fn(i,j')£m (*',/)) =aioc^exp(-X\n-m\)S iti fS jtj ', 

(4) 

respectively. Here a is the noise intensity of the global noise 
e n , &ioc is the noise intensity of the local noise £ n (hj), and 
A -1 is the correlation time of the local Gaussian colored noise. 
Here we set ui oc = a and A = 0.05. The parameter R mea- 
sures the noise correlation between a pair of neurons. r) n (i, j) 
is renewed at each iteration step n and for each unit in 
the iterated processing according to the algorithm proposed in 
Ref. |28| (see Appendix for details). In the following discus- 
sions, we will take a and R as controlled parameters. 



II. EQUATIONS OF THE NETWORK 

The Rulkov map2£i2! is employed to model the dynamical 
behavior of neurons constituting the examined neuronal net- 
work. The model captures succinctly main dynamical mecha- 
nisms in real neuronal ensembles, foremost showing typical 
restructuring of collective behavior following stochastic in- 
puts. Specifically, we consider a network of N x N electrically 
coupled Rulkov maps 

u n +i(i,j) = ot/[l +ul(i,j)] +v n (i,j) + D[u n (i + 1, j) 
+ u„(i - + u n (i,j - 1) + u n (i,j + 1) 
-4u„(i,j)] +T] n (i,j), 

Vn+l(i,j) = V n (i,j) ~ 0U n (i,j) - 7, 



HI. MEASURES OF NEURONAL DYNAMICS 

Two measures for quantifying the observed neuronal dy- 
namics due to the impact of noise are employed. One is the 
mean firing rate n, 29 and the other is the population coherence 
K i0 j 3i M Yjjg mean fi r i n g ra te of the neuronal network is defined 

as 

n = (ir(n)) T = /±^0[u n (iJ)-u t h]\ , (5) 

where Uth = —0.2 is the firing threshold determined by the 
action potential of the Rulkov neuron. Notably, 6(x) is a heav- 
iside function with 0(x) = 1 if x > and 9 = if x < 0. The 
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bracket () indicates the average over the whole iteration time 
T. 

To quantify the correlations of firing events in the neu- 
ronal network, we introduce the population coherence mea- 
sure k(t) i 30,31 k(t) is defined as the average of the local co- 
herence Kij (t) over all the pairs of neurons, namely: 



1 N 



N 



N(N-l) 

2 i,j=l\i^j 



where 



Kij(T) 



1=1 



(6) 



(7) 



i=i 



The coherence Kij(r) between any two neurons i and j is 
measured by the cross-correlation of their spike trains at zero 
time lag within a time interval t. More specifically, we divide 
the full iteration time T into small bins of duration t = 70 and 
define the two spike trains as Yi(l) = or 1 and Yj(l) = or 
1(1 = I, 2, ...,m;T/m = r), whereby Y(l) = 1 if the onset 
of a spike occurred at the Z-th time bin, otherwise Y(l) = 0. 
Sometimes, the population coherence measure k(t) is also 
used to quantify the synchronization of neuronal firings in 
networks.^Si^ Larger population coherence n corresponds to 
higher correlations between neurons inside the network. 



IV. EFFECTS OF CORRELATED GAUSSIAN NOISE 

In order to discern clearly the distinct effects of correlated 
Gaussian noise on neuronal dynamics, we consider firstly the 
effects of local Gaussian noise, i.e., R = 0.0. Variations of 
the mean firing rate II and the population coherence measure 
k with respect to the noise intensity a for various coupling 
strengthes D are shown in Figs. [JJ and [2] Importantly, local 
Gaussian noise is white in Figs.[JJbut colored in Figs. [2] From 
these two figures, we can see that II and k can reach larger 
values at an intermediate noise intensity than at smaller and 
larger noise intensities, for both local Gaussian white and col- 
ored noise. We caution, however, that for large noise intensi- 
ties the mean firing rate II may yield spurious results because 
of the very noisy output of individual neurons forming the 
neuronal network, due to which it is practically impossible to 
discern what is a firing event and what not. In accordance 
with this, the results depicted in this paper are constrained to 
noise intensities for which the neuronal dynamics still plays a 
significant role, i.e., is not completely overshadowed by noise. 

The above-reported results can be interpreted as follows. 
Small noise intensities are unable to evoke excitations, accord- 
ingly, the mean firing rate II becomes zero and the coherence 
of neurons inside the network is small (i.e. k is small). For 
intermediate noise intensities, typically only a few neurons (at 
random) forming the lattice start firing. Due to the diffusive 



coupling and the noisy support the excitations can propagate 
regularly to the neighbors, which results in well ordered circu- 
lar waves (as shown in Ref. [l3l) that ultimately result in high 
firing rate II and large coherence k. While for large noise in- 
tensities, neurons inside the network exhibit high-rate spiking 
behavior. When such spiking neurons are coupled diffusively, 
they tend to suppress the inputs coming from other neurons, 
which in turn decreases the number of firing events in a time 
span and destroys the spatial order of the dynamics, finally 
leading to small II and a decreases of k. Therefore, for local 
Gaussian noise, the mean firing rate and the correlations of 
the neuronal network can be optimized by some intermediate 
noise intensities, as shown in FigsQ]and[2] 

In what follows, we examine more closely the effects of 
correlated Gaussian noise via controlling the noise correlation 
R. Variations of the mean firing rate II and the population 
coherence measure k with respect to the noise intensity a for 
various noise correlations R are shown in Figsf3]and[4] The 
local Gaussian noise is taken as white in Figsf3] and colored 
in Figs0J respectively. Compared with the results shown in 
FigsQJa) and|2ja), we can see that variations of II with respect 
to the noise intensity a under correlated Gaussian noise are 
similar as the ones under local Gaussian noise. 

For the population coherence k, however, we can see that its 
variations versus the noise intensity a are strongly dependent 
on the noise correlation R as shown in Figsf3jb) and|4jb). It 
is also worth pointing out that the occasional non-smoothness 
of the curves in the later two figures is predominantly a conse- 
quences of the somewhat erratic switching between the emer- 
gence of spatially ordered patterns and their absence on the 
network. The latter introduces some non-smoothness to the 
employed statistical quantifiers, yet we found it impossible to 
eliminate this by means of more intensive numerical investi- 
gations. The reader is also referred to Ref. LL3J where pattern 
formation due to correlated Gaussian noise has been studied 
earlier. Through comparisons with the corresponding results 
presented in Figs. [Tib) and|2b), we find that effects of cor- 
related Gaussian noise on the population coherence k look 
more complex than the ones of local Gaussian noise. In the 
case of local Gaussian white noise, as shown in Fig. Ob), a 
weak coherent behavior (there exists a small peak at an inter- 
mediate a) can be observed for small R, e.g., R = 0.1,0.3. 
And for an intermediate R(R = 0.5, 0.7), the variations of k 
become irregular with changing of the noise intensity. With 
further increasing of R, a plateau region of n emerges as the 
noise intensity a increases, as can be observed from the line 
corresponding to R = 0.9, for example. For local Gaussian 
colored noise, variations of k with respect to a are irregular at 
R = 0.1, and exhibit plateau regions with increasing of a for 
large R, as presented in Fig.Hfb). 

Notably, the underlining mechanism regarding the effects 
of correlated Gaussian noise on the mean firing rate II is sim- 
ilar as discussed above in the case of local Gaussian noise. 
However, in order to clarify the effects of noise correlation R 
on the complicated variations of k versus er, we introduce the 
quantity 



2Rcr 



(8) 
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FIG. 1: Stack lines by Y offsets of the mean firing rate n (a) and the population coherence measure k (b) in dependence on the noise intensity 
a of additive local (we set R = 0.0) Gaussian white noise for various coupling strength D. Note that the 2>axis has a logarithmic scale. 




which is the noise strength ratio between global and local 
noise, p = R/(l — R) when local noise is taken as Gaus- 
sian white noise, and p = 2R/ (1 — R)X when it is taken as 
Gaussian colored noise. Thus, values of p can be controlled 
by the noise correlation R and A. We calculate the variations 
of k with respect to a for different noise intensity ratios p, or 
equivalently, different pairs of (R, A). The obtained results 
are presented in Fig. [31a), where pairs of R and A are colored 
black if the variation of k with respect to a shows a weak 
coherence. At this point it is instructive to examine the cor- 
responding curves for R = 0.1 and 0.3 that are depicted in 
Fig. 0b). Conversely, pairs of R and A are colored gray if 
variations of k with respect to a are irregular, as can be ob- 
served observed from the corresponding curves for R = 0.5 
and 0.7 in Fig. [2b) and from the curve depicted for R = 0.1 
in Fig. Htb). Finally, (R, A) pairs are colored white if there 
exists a plateau region of k for large a, as can be observed in 



Fig. 0b) for R = 0.9 and in Fig. Bib) for R = 0.3, 0.5, 0.7 
and 0.9. In order to now appreciate the quantity p introduced 
in Eq. (8) as an important driving force behind the variations 
k with respect to a, we show in Fig. |5jb) how p varies for 
different pairs of R and A in a systematic manner. In partic- 
ular, from thus far presented results it can be concluded that 
there exist three intervals of p, i.e. (0, a), [a, b), [b, oo), within 
which the variations of k with respect to a are different from 
one another. In Fig. |5Jb) black is used for those combina- 
tions of R and A for which 0.0 < p < 2.0, gray is used if 
2.0 < p < 8.0, while white is used if 8.0 < p < oo. Com- 
pared to results presented in Fig. EJa), we can observe at a 
glance that by setting a and b to equal 2.0 and 8.0, respec- 
tively, the color patterns match nearly perfectly, from which 
we conclude that p indeed has a decisive impact on variations 
of k with respect to a. 

From the analysis of the effects of correlated Gaussian 
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FIG. 3: Stack lines by Y offsets of the mean firing rate II (a) and the population coherence measure k (b) in dependence on the noise intensity 
a of additive correlated global noise r\(i, j) with local Gaussian white noise for various noise correlations 7?. The coupling strength is constant, 
equalling D — 0.0025. Note that the 2>axis has a logarithmic scale. 




FIG. 4: Stack lines by Y offsets of the mean firing rate n (a) and the population coherence measure /t (b) in dependence on the noise intensity a 
of additive correlated global noise rj(i,j) with local Gaussian colored noise for various noise correlations R. The coupling strength is constant, 
equalling D — 0.0025. Note that the a;-axis has a logarithmic scale. 



noise, we now thus know that noise correlations have no no- 
table effects on the variations of II versus a, while conversely, 
the noise correlation R plays a crucial role in how k varies 
with respect to a, in particular by means of controlling the 
noise strength ratio p. 



V. SUMMARY 

In this paper, we have studied effects of correlated Gaussian 
noise on neuronal firings, measured by the mean firing rate 
and the population coherence, of a two-dimensional network, 
which is locally modeled by the Rulkov map. Based on our 
numerical simulations, we have found that the mean firing rate 
of the network II can be enhanced at some intermediate noise 



intensities by correlated Gaussian noise for any noise corre- 
lation R. This phenomenon is similar as the coherence reso- 
nance, even though there is not a well-defined optimal noise 
intensity. While for population coherence measure k, we find 
that its variations with respect to noise intensity are very com- 
plex. Furthermore, we reveal that it strongly depends on the 
ratio p. Moreover, through the measure for population coher- 
ence k, we have been able to gain a deeper understanding of 
the interactions between global and local noise. In particular, 
we have shown that an appropriately tuned global noise can 
be an effective promotor of correlations of firing events in the 
neuronal network. 

As we have already stated in the Introduction, it is still de- 
batable whether a cortical neuron is driven mainly by the mean 
firing rate of presynaptic neurons or by correlations between 
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FIG. 5: (a) Variations of k with respect to a for different (R, A) pairs. If the coherence of k with respect to a is weak then the corresponding 
(R, A) pairs are colored black. If the variations of k with respect to a are irregular the color is gray, while it is white if there exists a plateau 
region of k for large a. (b) The dependence of p on different (R, A) pairs. Color black is used for those combinations of R and A for which 
0.0 < p < 2.0, gray is used if 2.0 < p < 8.0, while white is used if 8.0 < p < oo. See also the main text for further details. 



pairs of neurons inside the neuronal network. Moreover, neu- 
rons inside neuronal networks are not only affected by local 
random fluctuations, but also stimulated by some common 
random inputs. Thus, the results obtained for discussing ef- 
fects of correlated Gaussian noise on the mean firing rate and 
correlations of the neuronal network may give some important 
implications on investigating transmission of neuronal infor- 
mation in neuronal networks. 
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Appendix: Noise generation 

Gaussian white noise e„ can be generated effectively for 
map-based neuronal networks as follows. Let 

a = random number, 

b = random number, (A.l) 
e n = [-4o-ln(a)] 1 / 2 cos(2vr&), 



where a and b are uniformly distributed on the unit interval. 
When £, n {i,j) is Gaussian white noise, e n will be renewed 
according to Eq. ( A. 1 ) at each iteration step n and for each unit 
(i, j). When £ n (i,j) is a Gaussian colored noise, however, it 
can be generated by means of 



a = random number, 

b — random number, 

g w = [-4crAiln(a)] 1 / 2 cos(2vr6), 

£n+l = £n — A£ n + Xg w , 



with initial conditions 



I = random number, 

m = random number, (A. 3) 

g c = [-4crAln(0] 1/2 cos(27r77i), 



where a, 6, 1 and m are uniformly distributed on the unit in- 
terval. Subsequently, rj n (i,j) — \f~Re n + — R£ n {i,j) is 
renewed at each iteration step n and for each unit by 
repeating Eqs. (A. 2) and (A. 3). 
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